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Abstract —In this letter, we propose a practical non-linear energy 
harvesting model and design a resource allocation algorithm for 
simultaneous wireless information and power transfer (SWIPT) 
systems. The algorithm design is formulated as a non-convex 
optimization problem for the maximization of the total harvested 
power at energy harvesting receivers subject to minimum re¬ 
quired signal-to-interference-plus-noise ratios (SINRs) at multiple 
information receivers. We transform the considered non-convex 
objective function from sum-of-ratios form into an equivalent 
objective function in subtractive form, which enables the derivation 
of an efficient iterative resource allocation algorithm. In each 
iteration, a rank-constrained semidelinite program (SDP) is solved 
optimally by SDP relaxation. Numerical results unveil a substantial 
performance gain that can be achieved if the resource allocation 
design is based on the proposed non-linear energy harvesting 
model instead of the traditional linear model. 

I. Introduction 

Energy harvesting (EH) is a promising solution for prolong¬ 
ing the lifetime of communication networks by introducing 
self-sustainability to energy-limited devices. Among different 
EH technologies, wireless power transfer (WPT) via electro¬ 
magnetic waves in radio frequency (RE) enables comparatively 
controllable EH at the receivers compared to conventional 
natural energy sources, such as wind and solar. Recent progress 
in the development of RE-EH circuitries has made RE-EH prac¬ 
tical for low-power consumption devices, e.g. wireless sensors. 
In particular, RE-EH provides the possibility of simultaneous 
wireless information and power transfer (SWIPT) ID-ia. Yet, 
this new technology introduces a paradigm shift in system 
and resource allocation algorithm design. In 13 , the authors 
studied rate-energy trade-off regions by designing an optimal 
beamformer. In a, energy-efficient SWIPT was investigated in 
multicarrier systems, where power allocation, user scheduling, 
and subcarrier allocation were considered. In a, the authors 
solved the energy efficiency maximization problem for large- 
scale multiple-antenna SWIPT systems. However, existing lit¬ 
erature studies Ul and resource allocation algorithm designs for 
SWIPT networks are based on a linear EH model where 

the RE-to-direct current (DC) power conversion efficiency is 
independent of the input power level of the EH circuit. In 
practice, EH circuits 0-111 usually result in a non-linear end- 
to-end wireless power transfer. Hence, the conventional linear 
EH model cannot properly model the power dependent EH 
efficiency which leads to a mismatch for resource allocation. To 
the best of the authors’ knowledge, a practical non-linear EH 
model and a corresponding resource allocation algorithm design 
for SWIPT networks has not been reported in the literature, yet. 

In this letter, we address the above issues. To this end, we 
first propose a practical parametric non-linear EH harvesting 
model and verify its accuracy with measurement data. Then, 
we formulate the resource allocation algorithm design as a non- 
convex optimization problem for maximization of the total har¬ 
vested energy. The considered non-convex optimization problem 
is solved optimally by an iterative algorithm. Simulation results 
illustrate the total harvested energy loss when a conventional 
linear EH model is adopted for resource allocation algorithm 
design. 

Notation: A^, Tr(A), A~^, and Rank(A) represent the 
Hermitian transpose, trace, inverse, and rank of matrix A, 
respectively; A ^ 0 indicates that A is a positive semidefinite 
matrix; matrix I^r denotes the N x N identity matrix. 
denotes the space of all A x M matrices with complex entries. 



Fig. 1. A downlink SWIPT system with K = 2 information receivers (IRs) 
and J = 2 energy harvesting receivers (ERs). 

represents the set of all N-hy-N complex Hermitian 
matrices. The distribution of a circularly symmetric complex 
Gaussian (CSCG) vector with mean vector x and covariance 
matrix S is denoted by CA/'(x, S), and ~ means “distributed 
as”. £{■} denotes statistical expectation. 

II. System Model 

A. Channel Model 

We focus on a frequency flat slow fading channel for 
downlink multiuser SWIPT systems, cf. Eigure [T] In particular, 
a transmitter equipped with At > 1 antennas serving K 
information receivers (IRs) and J energy harvesting receivers 
(ERs) is considered. The K IRs are low complexity single¬ 
antenna devices and each ER is equipped with Ar receive 
antennas to facilitate EH. In each time slot, the transmitter sends 
a vector of data symbols to the K IRs. The received signals at 
IR k and ER j are given by 

K 

Vk = '^^kSk + nk, yke{l,...,K}, and (1) 

k^l 

K 

YeRj = Gf y^ WkSk + riERj, Vj g {1,..., J}, (2) 

k=l 

respectively, where Sfc S C and S are the 

data symbol and the beamforming vector intended for IR 
k, respectively. Without loss of generality, we assume that 
= IjVfc G A}. The channel vector between 

the transmitter and IR k is denoted by hj; G and the 

channel matrix between the transmitter and ER j is denoted by 
G^. e _ CA/'(0 ,(t 2) and iier, - CAr(0, 

are the additive white Gaussian noises (AWGN) at the IRs and 
the ERs, respectively, denotes the noise power at the receiver. 

B. Energy Harvesting Model 

In the literature, the total harvested energy at ER j, 
is typically modelled by the following linear model 13 - 0 : 

K 

‘^ERf = ^^^er,, Per, = ^Tr (w^wf G.Cf), (3) 

where PeRj is the received RF power at ER j and 0 < < 1 

is the fixed energy conversion efficiency of ER j. We note that 
in this linear EH model, the energy conversion efficiency is 
independent of the input power level at the ER. In other words, 
the total harvested energy at the ER is linearly and directly 
proportional to the received RE power. However, in practice, 
EH circuits 0-0 result in a non-linear end-to-end wireless 
power transfer. 





In general, for low power, the RF energy conversion effi¬ 
ciency improves as the input power rises, but there are diminish¬ 
ing returns and limitations on the maximum possible harvested 
energy, as was evidently proved by field measurements Ei-il. 
Thus, it is expected that the conventional linear EH model is 
only accurate for the specific scenario when the received powers 
at all ERs are constant. In this letter, we propose a practical 
parametric non-linear EH model which captures the dynamics 
of the RE energy conversion efficiency for different input power 
levels. In order to isolate the system model from the specific im¬ 
plementation details of the EH circuit, we propose a non-linear 
EH model based on the logistic (sigmoidal) function. Hence, 
the total harvested energy at ER j, modelled as; 
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, (4) 


(5) 


Here, vj/P^^ticai traditional logistic function with respect 

to the received RF power Per^ - We introduce a constant 
in (HI to ensure a zero-input/zero-output response for EH. Mj 
is a constant denoting the maximum harvested power at ER 
j when the EH circuit is saturated. Parameters aj and bj are 
constants related to the detailed circuit specifications such as the 
resistance, capacitance, and diode turn-on voltage. In practice, 
the EH hardware circuit of each ER is fixed and the parameters 
Qj, bj, and Mj of the proposed model in (|4]i can be easily found 
by a standard curve fitting tool. We note that the proposed 
non-linear EH model is able to capture the joint effect of the 
non-linear phenomena caused by hardware constraints including 
circuit sensitivity limitations and current leakage Q, ID. 

Eigure |2] illustrates that the proposed non-linear EH model 
closely matches experimental results 111, HI for the wireless 
power harvested by practical EH circuits. Besides, Eigure|2]also 
shows that the linear model in Q is not accurate in modelling 
non-linear EH circuits. 

In the sequel, we adopt the proposed non-linear EH model 
for resource allocation algorithm design. We assume that perfect 
channel state information is available for resource allocatiorO. 
Eurthermore, since does not affect the design of beam- 
formii^ vector w^, cf. (0, for simplicity, we will directly 
use jQ represent the harvested power at ER j in the 

following study. 


HI. Problem Eormulation and Solution 

The system design objective is the maximization of the total 
harvested power and can be mathematically formulated as: 


maximize 

Wfc 


subject to 


E ijfPractical 

^er, 

i=i 

K 

Cl : ^llWfclP < Pmax, 


(6) 




(^2- w"HfcWfe 


wf HfeWj- -b al 




where Hfc = hfeh^. Constants Pmax and F™'^ in constraints 
Cl and C2 are the maximum transmit power for the trans¬ 
mitter and the minimum required signal-to-interference-plus- 
noise ratio (SINR) at IR k, respectively. It can be observed 
that the objective function in ® is in the form of sum-of- 
ratios which is a non-convex function. In the following, we 
assume that the considered optimization problem is feasible for 


Fig. 2. A comparison between the harvested power for the proposed model 
in E) and the measurement data from two different practical EH circuits with 
different dynamic ranges from (7) and (8). The parameters aj, bj, and Mj in 
E) are calculated by a standard curve fitting tool. 

TABLE I 

Iterative Resource Allocation Algorithm. 


Algorithm Iterative Resource Allocation Algorithm 

1: Initialize the maximum number of iterations Lmax, iteration index n = 0, 
fi, and /3 

2: repeat {Outer Loop} 

3: Solve the inner loop problem in ED via SDP relaxation for given 

obtain the intermediate beamformer 
4: if fl3l is satisfied then 

5: return Optimal beamformer wf = w', 

6: else 

7: Update and /3 according to ED and n = n + 1 

8: end if 

9: until 1131 is satisfied or n = Lmax 


the study of resource allocation algorithm design. Although the 
Dinkelbach method or the Charnes-Cooper transformation 
can be exploited to handle a single-ratio objective function, 
they cannot be applied for a sum-of-ratios objective function. In 
order to obtain a tractable solution, we first transform the non- 
convex objective function into an equivalent objective functiorQ 
in subtractive form via the following theorem. 

Theorem 1: Suppose is the optimal solution to 
then there exist two vectors fi* = [p*,..., p j] and /3* = 
[PI,, P*j\ such that is an optimal solution to the following 
optimization problem 


J 

maximize 


wf G.?- 


/ ^ I 


Mj-P* (l -b exp (-ai(-PER, -6,))) 


(7) 


where T is the feasible solution set of ®. Besides, also 
satisfies the following system of equations: 

/3;(l + exp(-a,(P^R^-6,))) -M, = 0, (8) 

/i*(l + exp(-a,(P,^R^.-6,))) -1 = 0, (9) 

and = Ef=i Tr (w^(w^)«G,Gf). 

Proof: Please refer to ifTOl . ifTTl for a proof of Theorem 1. 

Theorem [T] suggests that for the maximization problem with 
sum-of-ratios objective function in ®, there exists an equivalent 
parametric optimization problem with an objective function in 
subtractive form, such that both problems have the same optimal 
solution w^. As a result, the optimization problem can be solved 
by an iterative algorithm consisting of two nested loops. In the 
inner loop, we solve the optimization in 0 for given (/x,/3). 
Then, in the outer loop, we find the optimal (/r*,/3*) satisfying 
the system of equations in ® and (|9]l, cf. algorithm in Table U 


* In practice, the channel state information of the IRs/ERs can be obtained 
during the handshaking between the transmitter and the receivers before power 
and data transfer start. 


^Here, “equivalent” means that the optimization problem with the trans¬ 
formed objective function leads to the same resource allocation policy as the 
original problem. 






















A. Solution of the Inner Loop Problem 

As shown in Table |I] in each iteration in the inner loop, i.e., 
in line 3, we solve the following optimization problem for given 
parameters {11,(3): 


maximize > u, M,- — /3,- (1 + exp f — a,-(r, — bi 

WfceH«T.r, ^ ^ 

j=i 

K 

subject to Cl : ET>-(Wfe) Cmax; 

fc=l 

Tr(HfeWfc) 


( 10 ) 


C2 : 


preq 


>ETr(HfcW,)+a2,Vfc. 


C3: Rank(Wfc) <1,VA:, 


K 

C4 : r, < E Tr (WfcG.Gf), Vj, C5 : ^ 0, Vfc, 

fc=i 


where and Tj are the new and auxiliary optimiza¬ 

tion variables, respectively. Although the transformed objective 
function is in subtractive form, the transformed optimization 
problem in (fTOl i is still non-convex due to the rank-one matrix 
constraint C3. To obtain a tractable problem formulation, we 
apply SDP relaxation. Specifically, we relax constraint C3 in 
Coll by removing constraint Rank(Wfc) < 1 from the problem. 
Then, the considered problem becomes a convex SDP problem 
and can be solved by standard numerical algorithms for convex 
programs such as the interior point method. Now, we study the 
tightness of the SDP relaxation in the following theorem. 

Theorem 2: Assuming that the channels, i.e., hfc and Gj, 
are statistically independent and (fTol) is feasible, the optimal 
beamforming matrix of the SDP relaxed problem of (fTol i is a 
rank-one matrix with probability one, i.e., Rank(Wt) = 
for > 0. 

Proof: Please refer to the Appendix. 

Therefore, the adopted SDP relaxation is tight whenever the 
general channel conditions stated in Theorem 2 are satisfied. 
Hence, beamforming is optimal for the maximization of total 
harvested power for the proposed non-linear EH model. 


B. Solution of the Outer Loop Problem 

In this section, we present an algorithm to iteratively update 
(fi, (3) for the outer loop problem. For notational simplicity, we 

define functions = /3j(^l + exp ( —aj(PERj — 

Mj and + exp(-ai(PERi - - 1, 

i G J}. It is shown in ifTOl . ifTTl that the unique 

optimal solution {ij,*,(3*) is obtained if and only if (p{fi,(3) = 
Vpi: Pii ■ ■ • I P 2 j] = 0- Thus, the well-known damped Newton 
method can be employed to update (fi, (3) iteratively. In partic¬ 
ular, in the n-th iteration, and /3”+^ can be updated as, 

respectively, 

= /.t”"-F C"q" and /3"+i =/S’"-F (11) 

where q” = /S)]" /3) (12) 

and <p'(/r, (3) is the Jacobian matrix of (3). Cf’ is the largest 
e* satisfying 

||(p(/r’"+e'q",/3"+eV)ll (13) 

where I G {1, 2,...}, e* G (0,1), and S G (0,1). The damped 
Newton method converges to the unique solution {ij,*,(3*) 
satisfying the system of equations ® and (|9l), cf. ifTOll . ifTTI . 


Remark 1: We note that when there is one ER in the system, 
the traditional linear and the proposed non-linear EH model will 
lead to the same optimal resource allocation policy. 

Remark 2: The signal model adopted in this paper can be 
extended to include dedicated energy beams for the ERs by 
following a similar approach as in Ha. 



Fig. 3. Average total harvested power (dBm) versus the minimum required 
SINR (dB). 



Fig. 4. Average total harvested power (dBm) versus the number of ERs. 


IV. Results 

In this section, we present simulation results to demonstrate 
the system performance of the proposed resource allocation 
algorithm design. We assume a carrier center frequency of 
915 MHz and a signal bandwidth of 200 kHz. There are 
K = 2 IRs and J ERs located 50 meters and 10 meters 
from the transmitter, respectively. Each ER is equipped with 
Vr = 2 receive antennas. Since the ERs are close to the 
transmitter, line-of-sight communication channels are expected. 
Hence, the multipath fading between the transmitter and the ERs 
is modeled as Rician fading with a Rician factor of 3 dB. In 
contrast, the IRs are located 50 meters from the transmitter, thus, 
a line-of-sight may not be available and the multipath fading 
between the transmitter and the IRs is modeled as Rayleigh 
fading. All transmit antenna gains are 10 dBi. The thermal noise 
power is = —95 dBm. Eor the non-linear EH circuits, we set 
Mj = 20 mW which corresponds to the maximum harvested 
power per ER. Besides, we adopt Uj = 6400 and bj = 0.003 
which were obtained by curve fitting for measurement data from 
Q. The average system performance is obtained by averaging 
over different channel realizations. 

Eigure [3] depicts the average total harvested power versus 
the minimum received SINR at the IRs for J = 10 ERs and 
different numbers of transmit antennas. We assume that all IRs 
require the same minimum receive SINR, i.e., T™'^ = T"®'! 
and Pinax = 30 dBm. Extensive simulations (not shown 
here) have revealed that, in general, the proposed iterative 
algorithm converges to the globally optimal solution after less 
than 10 iterations. It can be observed from Figure [3 that the 
average total harvested power is a monotonically decreasing 
function with respect to Indeed, to satisfy a more stringent 
minimum SINR requirement, the transmitter is forced to steer 
the direction of transmission towards the IRs leading to a 
smaller amount of RF energy for EH. On the other hand, 
the total harvested energy increases for an increasing number 
of transmit antennas Vt, since the extra degrees of freedom 
offered by the increased number of transmit antennas facilitates 



















a more power efficient resource allocation. For comparison, we 
also show the performance of a baseline scheme in Figure [3 
For the baseline scheme, the resource allocation algorithm is 
optimized for maximization of the total system harvested power 
according to the conventional linear EH model in Q subject to 
constraints Cl and C2. As can be observed, the baseline scheme 
can only achieve a strictly smaller amount of total harvested 
power due to the resource allocation mismatch. In particular, 
the baseline scheme may cause saturation in EH in some ERs 
and underutilization of other ERs because it does not account 
for the non-linear nature of the EH circuits. 

Eigure |4] shows the average total harvested power versus the 
number of ERs J for Pmax = 30 dBm, a minimum required 
SINR of 30 dB, and different numbers of transmit antennas 
Nt. It can be observed that the average total harvested power 
increases with the number of ERs and the number of transmit 
antennas. In fact, a larger portion of the radiated power can 
be harvested when there are more ERs in the system since 
more receivers participate in the EH process. Besides, the 
performance gain of the proposed scheme compared to the 
baseline scheme increases with increasing number of ERs. This 
is because the resource allocation mismatch for the baseline 
scheme becomes more pronounced for a larger number of ERs 
leading to unsatisfactory performance. 

V. Conclusions 

In this letter, we proposed a practical EH model to capture 
the non-linear characteristics of EH circuits in SWIFT systems. 
Eurthermore, the resource allocation algorithm design for the 
proposed model was formulated as a non-convex optimization 
problem with a sum-of-ratios objective function and was solved 
optimally by the proposed iterative algorithm. Our simulation 
results unveiled that resource allocation algorithms designed 
for the conventional linear EH model, which is widely used 
in the literature, may lead to resource allocation mismatches 
for practical non-linear EH circuits. 

Appendix-Proof of Theorem[2] 

It can be verified that strong duality holds for the SDP relaxed 
version of (doll. Thus, solving the dual problem is equivalent to 
solving the primal problem da. Now, we prove Theorem |2] by 
first defining the Lagrangian function; 

K K TT 

C =^Tr(BfcWfe) +X! Tr ((Y^ + + A, (14) 

J 

Bfc =-AIjv, - ^ 7,H, + E ftG.Gf ^ (15) 

j=i 

where A is the collection of all constants and variables that 
are independent of and are thus not relevant in the 

proof. A, 7 fc,/ 9 j, and Y^ are the dual variables associated with 
constraints Cl, C2, C4, and C5, respectively. Then, the dual 
problem of the SDP relaxed problem of (fTOl i is given by 

minimize maximize C. (16) 

A,7fc.pj>0,Yfc^O 

We define and H* = {A*, 7 *,pj, Y*} as the set of 

optimal primal and dual variables of the SDP relaxed version 
of ( [Tol l, respectively. Now, we consider the following Karush- 
Kuhn-Tucker (KKT) conditions; 

Y*^0, A*, 7 :,P* >0,Vfc, Vj, (17) 

Y^W^ = 0 , ^ (18) 

Y^ = -B^-^, (19) 

where B*. is obtained by substituting the optimal dual variables 
H* into (Il5b . Equation ( fTSl l is the complementary slackness 
condition which is obtained by taking the derivative of the 


Lagrangian function with respect to Y^. Besides, equation ( fTSl l 
indicates that the columns of lie in the null space of Y^ 
for 7 ^ 0. Therefore, if Rank(Y^) = — 1, then the 

optimal beamforming matrix is a rank-one matrix. To reveal 
the structure of Y^, we show by contradiction that BJ is a 
negative definite matrix with probability one. Eor a given set of 
optimal dual variables, H*, and r*, (fTfil) can be written as 

maximize C. (20) 

Wfc gH~‘t 

Suppose B^ is not negative definite, then we can construct 
Wfe = rVfeV^ as one of the optimal solutions of (l 20 l i. where 
r > 0 is a scaling parameter and is the eigenvector 
corresponding to one of the non-negative eigenvalues of B^. 
We substitute into ( l20l ) which leads to £ = 

T,k=i Tr(LB* vfevf) -f r J2k=i Tr (vfevf (Y* -f ^)) + A. 
Since the channels of Gj and hfc are assumed to be statistically 
independent, it can be shown that 7 ^ > 0 for the optimal 
solution. Also, it follows that by setting r ^ 00 , the dual 
optimal value becomes unbounded from above. However, the 
optimal value of the primal problem is finite for a finite 
Pmax- Thus, strong duality does not hold which leads to a 
contradiction. Therefore, B^ is a negative definite matrix with 
probability one, i.e., Rank(B^) = At. Then, by exploiting (fT^ 
and a basic inequality for the rank of matrices, we have 

Rank(-Y^) = Rank(Y^) = Rank (b^ + 

, *TT , 

> Rank(B^) - Rank = At - 1. (21) 

Eurthermore, 7 ^ 0 is required to satisfy the minimum SINR 
requirement of IR k in C2 for F^®'^ > 0. Hence, Rank(Y^) = 
At - 1 and Rank(W^) = 1. ■ 
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